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It is shown that the KadanofT-Baym equations at consistent first-order 
gradient approximation reveal exact rather than approximate conservation 
laws related to global symmetries of the system. The conserved currents and 
energy-momentum tensor coincide with corresponding Noether quantities in 
the local approximation. These exact conservations are valid, provided a <1> 
derivable approximation is used to describe the system, and possible memory 
effects in the collision term are also consistently evaluated up to first-order 
gradients. 



1. INTRODUCTION 

Non-equilibrium Green function techniques, developed by Schwinger, 
KadanofF, Baym and Keldysh [1-4], provide the appropriate concepts to 
study the space-time evolution of many-particle quantum systems. This 
formalism finds now applications in various fields, such as quantum chromo- 
dynamics [5,6], nuclear physics, in particular heavy ion collisions [7-23], as- 
trophysics [11,24,25], cosmology [26], spin systems [27], lasers [28], physics 
of plasma [29,30], physics of liquid '^He [31], critical phenomena, quenched 
random systems and disordered systems [32], normal metals and super- 
conductors [24,33,34], semiconductors [35], tunneling and secondary emis- 
sion [36], etc. 

For actual calculations certain approximation steps are necessary. In 
many cases perturbative approaches are insufficient, like for systems with 
strong couplings as treated in nuclear physics. In such cases, one has to 
resum certain sub-series of diagrams in order to obtain a reasonable approx- 
imation scheme. In contrast to perturbation theory, for such resummations 
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one frequently encounters the fact that the scheme may no longer be con- 
serving, although for each diagram considered the conservation laws are 
implemented at each vertex. Thus, the resulting equations of motion may 
no longer comply with the conservation laws, e.g., of currents, energy and 
momentum. This is a problem of particular importance for recent stud- 
ies of particles with broad damping width such as resonances [18]. The 
problem of conservation laws in such resummation schemes has first been 
considered in two pioneering papers by Baym and Kadanoff [37,38] dis- 
cussing the response to an external perturbation of quantum systems in 
thermodynamic equilibrium. Baym, in particular, showed [38] that any 
approximation, in order to be conserving, must be based on a generat- 
ing functional $. This functional was first considered by Luttingcr and 
Ward [39] in the context of the thermodynamic potential, cf. [40] , and later 
reformulated in terms of path integrals [41]. For truncated self-consistent 
Dyson rcsummations this functional method provides conserved Noether 
currents and the conservation of total energy and momentum at the ex- 
pectation value level. In our previous paper [15] we extended the concept 
to the real-time Green function technique and relativistic systems, con- 
structing conserved 4-currents and local energy-momentum tensor for any 
chosen approximation to the $ functional. While ^-derivable Dyson re- 
summations formulated in terms of the integro-differential Kadanoff-Baym 
(KB) equations indeed provide exact conservation laws, these equations 
are usually not directly solvable. Therefore many applications involve fur- 
ther approximations to the KB equations: the gradient approximation and 
often the quasi-particle approximation leading to differential equations of 
mean field and transport type. Any improvement of the quasi-particle ap- 
proximation beyond the mean-field level, e.g., through inclusion of energy- 
momentum-dependent self- energies, again leads to difficulties with conser- 
vation laws. Various attempts to remedy this problem were undertaken, 
see refs. [35,42-47] and references therein. An essential progress within the 
quasi-particle approximation was achieved by the ansatze of refs. [35,44]. 

Interested in the dynamics of particles with broad mass width like res- 
onances we like to discuss the question of conservation laws for transport 
problems at a much more general level, than usually considered. We call 
this the quantum transport level. It completely avoids the quasi-particle 
approximation, and solely rests on the first-order gradient approximation 
of the KB equations. This concept was first addressed by Kadanoff and 
Baym [2] in the chapter "Slowly varying disturbances", Eqs. (9-25). Moti- 
vated by applications for the description of heavy-ion collisions further at- 
tempts were recently suggested [18-23] based on the so called Botermans 
Malfliet (BM) substitution [10]. Within the BM choice of the quantum 
kinetic equations a number of desired properties including an H-theorem 
for a local entropy current related to these equations were derived [18], 
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however no strict realization of the conservation laws for the Noether cur- 
rents were obtained [18,20]. Due to the approximation steps involved one 
may expect the quantum transport equations both, in the KB and BM 
forms, to possess only approximate conservation laws though in line with 
the level of approximation. Such approximate nature of conservation laws 
may be well acceptable theoretically. Nevertheless, both from a principle 
perspective and also from a practical point of view this situation is less 
satisfactory. Quantum kinetic equations which possess exact conservation 
laws related to the symmetries of the problem could serve as a natural ex- 
tension of the quasiparticle transport phenomenology to broad resonances, 
e.g. applicable to high energy heavy ion collisions. 

In this paper we give a proof that the quantum kinetic equations in the 
form originally derived by Kadanoff and Baym in fact possess the generic 
feature of exact conservation laws at the expectation value level. This holds 
provided all self-consistent self-energies are generated from a 4>-functional 
and all possible memory effects due to internal vertices within the self- 
energy diagrams are also consistently expanded to first-order gradients. 

In sect. 2 we review the derivation of the quantum kinetic equations 
in a notation suitable for our later derivation of the conservation laws in 
sect. 3. Sect. 4 deals with the general gradient approximation and its 
representation in terms of diagrams. Finally, in sect. 5 within the 
derivable method we give a diagrammatic proof of the exact conservation 
laws of the quantum kinetic equations for the KB choice. 

We restrict the presentation to physical systems described by complex 
quantum fields of different constituents interacting via local couplings. The 
kinematics can be either relativistic or non-relativistic. Extension to real 
boson fields, as well as to relativistic fermions is straight forward though te- 
dious in the latter case. We also exclude theories with derivative couplings. 



2. KADANOFF-BAYM EQUATIONS AND COMPLETE 
GRADIENT APPROXIMATION 

We assume the reader to be familiar with the real-time formulation of 
non-equilibrium field theory on the so called closed time contour, Fig. 1. 
Since we will deal with general multi-point functions we use the more con- 
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FIG. 1. Closed real-time contour with two external points x, y. 
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venient { — h} contour- vertex notation of refs. [4, 18]. Some rules are sum- 
marized in Appendix A.l. The set of coupled KB equations on the time 
contour in — h notation^ reads 

(G-i(-iai) - G^H-id2)) G-+il, 2) 

d3 (5](1-, 3)G(3, 2+) - G(l-, 3)S(3, 2+)) = tC(1-, 2+) (1) 



ic 

with Fourier transform of the inverse free Green function 

Q-i/ \ { ^ n^"^ for relativistic bosons 

° \po~"P^/(2to) for non-rel. fermions or bosons. ^ ' 

Here and below the upper/lower signs refer to fermions or bosons respec- 
tively, Go and G correspondingly denote the free and full Green functions, 
labels for the different species and internal quantum numbers are sup- 
pressed. The driving term on the r.h.s. of Eq. (1), summarized by C, 
is a functional of the Green functions through the self-energies E contour 
folded with G. The real-time integration contour, cf. Fig. 1, is denoted 
by C. The step towards transport equations is provided by introducing the 
four- dimensional Wigner transforms for all two-point functions through 

n^^y) = J (3) 
The KB Eq. (1) then transforms to 

v'^d^{Ti)G-+{X,p) = C-+{X,p;{G}) with t;/^ = Ag-1(p), (4) 

where the r.h.s. is also expressed in terms of the Wigner transforms of all 
Green functions through (3). The final step is to expand the complicated 
r.h.s. of Eq. (4) to the first-order gradients. Then the local part of this 
r.h.s., G^Yo^-j, consists of non-gradient terms, where one replaces the different 
mean positions {xi + Xj)/2 occurring in the various Green functions by the 
externally given mean position X of the l.h.s., i.e. X = + X2)/2 and 
evaluates the diagrams as in momentum representation. The corrections 



-"^The numbers 1, 2 and 3 provide short hand notation for spaee-time eoordinates 
XI, X2 and 2:3, respeetively, including internal quantum numbers. With superscript, like 
1~ and 2+, assigned to them, they denote contour coordinates with — and + specifying 
the placement on the time or anti-time ordered branch. Decomposed to the two branches 
the contour functions are denoted asF'='(l,2) = ^(1^,2') withfc,Z S {-,-|-}. The match 
to the notation used, e.g., in refs. [2,7] is given by F ^ = F^; F'^ = F-*; F = 
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for the displacement to the true coordinates of each Green function are then 
accounted for to the first order in the gradients. Here we simply abbreviate 
the gradient terms by a O operator acting on the local diagram expression 

v'^d^{Ti)G"+iX,p) = (l + iO){C(7+)(X,p)}, (5) 

where 

C(7+)(X,p) = C-+(X,p;{G(,o,)}), (6) 

is a functional of the local Green functions G(ioc) = G{X,p). Here and 
for all further considerations below, both, Green functions G{X,p) and 
sclf-cncrgics whenever quoted in their Wigncr function form, are 

taken in local approximation, i.e. with X given by the external coordinate 
and Y,{X,p) void of any gradient correction terms. The explicit definition 
of the O-operator is deferred to sect. 4. AH what we need to know at this 
level is that it consists of terms where in the diagrams defining C(ioc) pairs 
of Green functions are replaced by their space-time and momentum deriva- 
tives, respectively, just leading to equations linear in space-time gradients. 
Naturally the result of the O-operation depends on the explicit form, i.e. 
diagrammatic structure, of the functional on which it operates. Therefore 
for the non-gradient term in Eq. (5), which defines the local collision term 

= T^-\X,p)akiG'+{X,p) - {T)G-\X,p)aH^'+{X,p) (7) 

diagram 

= TiS-+(X,p)iG+-(X,p) - (T)iG-+(X,p)iI]+-(X,p) , (8) 

value ^ V ^ V 

gain loss 

we give both, the diagram expression (7) and the normally quoted value 
expression (8). The latter simplifies due to a cancellation of terms which 
however survive for the order sensitive gradient operation. Here cjik = 
a'^^ = diag(l, — 1) defines the "contour metric", which accoimts for the 
integration sense, and summation over the contour labels k,l € {— ,+} is 
implied, cf. (A.2) ff. 

The above quantum kinetic equation (5) has to be supplemented by a 
local Dyson equation for the retarded Green function [2] 

= (G^ip))-' - (9) 

which together with Eq. (5) provides the simultaneous solution to G"' . 
The full retarded Green function G^ depends on the retarded self-energy 
= T, — E = E+~ — E++ again in local approximation. Gq is 
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the free retarded Green function. Please note that equation (9) is just 
algebraic although it is obtained in the framework of the first-order gradient 

approximation. 

In most presentations of the gradient approximation to the KB equations, 
Eq. (5) is rewritten such that the gradient terms are subdivided into two 
parts, consisting of Poisson bracket terms describing drag- and back-flow 
effects, on the one side, and a memory collision term C'™''™, on the other 
side, in cases when the self-energy contains internal vertices 

v^'d^{T)iG-+{X,p) = {ReE^,TiG-+} + {TiS"+,ReG^} 

+G("nim)(^.rt + C(T+)(X,p), (10) 

where 

= S^^%(Xp)aHG'+(Xp)-G-'=(X,rta,,E|+^^)(X,p) (11) 

= -Srm'exn)(^>P)G+-(^,P)+G-+(X,p)S+^-^)(X,p) (12) 
value 

and 

J:imeu.){X,p) = ^OmX,p)}. (13) 

One may further introduce the spectral function A{X.p) = — 2ImG^(X,p) 
determined by the retarded equation (9), as well as the four-phase-space 
distribution function f{X,p), by means of =FiG ^{X.p) — f{X,p)A{X,p), 
whose evolution is governed by transport Eq. (10) or equivalently by (5). 
This defines a generalized quantum transport scheme which is void of the 
usual quasi-particle assumption. The time evolution is completely deter- 
mined by initial instantaneous values of the Green functions and their gra- 
dients at each space-time point, which means it is Markovian, since the 
memory part of the collision term is kept only up to first-order gradient 
terms. Within its validity range this transport scheme is capable to describe 
slow space-time evolutions of particles with broad damping width, such as 
resonances, within a transport dynamics, now necessarily formulated in the 
four-dimcnsional phase-space. 

Depending on the questions raised, the above separation (10) may not 
be always useful. For the derivation of the conservation laws only a uni- 
fied treatment of both, the Poisson brackets and memory collision terms, 
reveals the symmetry among these terms, which then displays the neces- 
sary cancellation of certain contributions such that the conservation laws 
emerge. Therefore, in the fortliconiing considerations we will mostly refer 
to the quantum kinetic equation in the form (5). 
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3. CONSERVATION LAWS 

Conservations of charge and energy-momentum result from taking the 
charge or four-momentum weighted traces of the transport equations (5). 

These traces inchidc the integration over four-momentum, as well as sums 
over internal quantum numbers and all types of species a with charges Ca 

The the charge and four-momentum leaks on the r.h.s, abbreviated as 

Q and T'^, can be represented by closed diagrams, where the two end 
points {x^ and X2 in Eq. (1)) coalesce, i.e. Xi = X2 = X. In coordinate 
representation the r.h.s. corresponds to contour integrals of the type (A. 7) 
- (A. 9) for which entirely retarded terms drop out. Therefore, one even 
can place the two end points on the same contour side (respecting the 
fixed order for Tad-pole terms). For definiteness we have chosen the time 
ordered (— ) branch. The external point x^ = .x^ = X is then the reference 
point with respect to which the gradients are to be evaluated. In line 
with causality requirements this reference point X = {t,x) is also the 
retarded point. This implies that in a real-time contour representation any 
contribution to (14) from internal integrations with physical times larger 
than t drops out, cf. [32,8]. Using the explicit form of the contour metric 
a (cf. Eqs. (A. 3) and (A. 4)) one obtains 

C;^ioc)i^^P) = T(S^fe(^,p)Gr(X,p)-G-'=(X,p)S^_(X,p)). (15) 

In the proof given in sect. 5 we show that the local parts of the r.h.s. of 
(14) as given by 

(tS)=?/lf.C">*.,..,.. 

entirely drop out. Also the gradient terms of Q cancel, while the gradients 
of the term compile to a complete divergence 

Q{X) = 0, T^(X) = g^'■'^|_, (£p°*(X) - £'^\X)) , (17) 

provided the self-energies are derived from a so called ^-functional [38]. 
This implies exact conservation laws for the Noether currents and the 
energy-momentum tensor given by 



d^j^'ix) = 0, a^er„';(x) = with 



(18) 
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J^{X) = Y,j ^,eaV^{Ti)G-+{X,p), (19) 

" +g^''{£\^l{X)-£f:^{X)). (20) 

Here G is the self-consistent propagator solving the coupled set of quantum 

transport equations (5) and (9). Furthermore, 0[^'^(X) is the Zoca/ version 
of energy-momentum tensor which for the ^-derivable approximation to 
the KB equations has been constructed in our previous papers [15,18]. 
Thereby, £™]:{X) and £y°^{X) define the interaction and single-particle 
potential energy densities, respectively, also taken in the local approxima- 
tion. 

In the subsequent sections we formally define the complete first-order 
gradient expansion for any two point function, which contains internal ver- 
tices, and specify the corresponding diagrammatic rules. Finally, in sect. 5 
wc prove the conservation laws (18), using the ^-derivable properties and 
the gradient rules. 



4. COMPLETE GRADIENT APPROXIMATION 

Let M(l,2) be any two-point function with complicated internal struc- 
ture. We are looking for its Wigner function M{X,p) with X = \{xi +X2) 
to first-order gradient approximation. The zero-order term is just given by 
evaluating M(l, 2) with the Wigner functions of all Green functions taken 
at the same space-time point X = (xi + X2)/'2. and the momentum inte- 
grations being done as in the momentum representation of a homogeneous 
system. To access the gradient terms related to any Green function G{i,j) 
involved in M(l,2), its Wigner function G{^{xi + Xj),p) is to be Taylor 
expanded with respect to the reference point X = {xi + X2)/2, i.e. 

G(^^,P) « G{X,p) + i [(xf - x'^) + {x^ - 4)] ^G{X,p). (21) 

Both the space derivatives of Green fimctions and the factors (xi —xi) and 
{xj — X2) accompanying them can be taken as special two-point functions, 
and we therefore assign them special diagrams 

= \{di + d^)G{i,j)^dxG{X,p), (22) 

— = -i(a;i-a;,) ^-(27r)4— (23) 

i j '^P 
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with the corresponding Wigner functions at the right hand side. Then 
the gradient terms of a compHcated two-point function (given in different 
notation) can graphically be represented by the following two diagrams on 
the r.h.s. 



0{M(1,2)} 



O 



= lOAf 
2 1 




3^ 



M' 



M' 



) (24) 



Here the diamond operator O, as above, formally defines the gradient 
approximation of the two-point function M to its right with respect to the 
two external points (1,2) displayed by full dots. The diagrammatic rules 
arc then the following. For any G(3, 4) in M, take the spatial derivative 
dxG{X,p) (double line) and construct the two diagrams, where external 
point 1 is linked to 3 by an oriented dashed line, and where point 2 is linked 
to 4, respectively. Interchange of these links provides the same resiilt. Here 
M' is a four-point function generated by opening M(l,2) with respect to 
any propagator G(3,4), i.e. 

^^^'^'^'^^ = ^^iG(4:^- ^ ^ 

The diagrams in Eq. (24) are then to be evaluated in the local approxi- 
mation, i.e. with all Wigner Green functions taken at the same space-time 
point X. The dashed line (23) adds a new loop integration to the diagram, 
which, if integrated, leads to momentum derivatives of the Green func- 
tions involved in that loop. Both, double and dashed lines have four-vector 
properties, and the rule implies a four-scalar product between them. 

The explicit properties of the dashed line (23) permit to decompose it 
into two or several dashed lines through the algebra 



1 



— — 1 — and • \ = 0. (26) 

1 2 2 3 ^ — 



In momentum-space representation these rules correspond to the partial 

integration. They imply that 0{M(1,2)} = 0, if M contains no internal 
vertices. Applying rule (26) to the convolution of two two-point functions 



C(l,2) = j d3A(l, 3)^(3,2) 



(27) 
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leads to the following convolution theorem for the gradient approximation 
0{C{X,p)}= o 

+ •^>^> + <^x;^ (28) 

= {A{X,p),B{X,p)} 

+A{X,p)0{B{X,p)} + 0{A{X,p)}B{X, P). (29) 

Besides the standard Poisson bracket expression {A, B} it leads to fur- 
ther gradients within each of the two functions (note that the O operator 
acts only on the two-point function in the immediate braces to its right). 
Applied to the r.h.s. of Eq. (5), this rule indeed provides the decompo- 
sition into Poisson bracket and memory terms for the more conventional 
formulation of the quantum kinetic equation (10). 

5. ^-DERIVABLE SCHEME AND EXACT CONSERVATION 

LAWS 

In this section we first review the ^-derivable properties at the level 
of self-consistent Dyson or KB equations, then proceed towards the im- 
plications for the quantum kinetic equations (5) and to the proof of the 
corresponding exact conservation laws (18). 

In a <i>-derivable scheme the self-energies are generated from a functional 
${G, A} through the following functional variation^, cf. [15] 

_ Si<^{G,\} ( 2 for real fields , , 

il.[x,y) - T^^^^— ^x|^ for complex fields " ^'^^> 

The $ functional itself is given by two-particle irreducible {'2.PI) closed 
diagrams in terms of full Green functions of the underlying field theory, i.e. 



\ \ Jc 



i${G, A} = ( exp i / d'*a:A(a;)£'"* ) . (31) 



2PI 



^Here we include also the rule for real fields, upper /lower signs refer to 
fermions /bosons 
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The interaction strength A(a;), the physical value of which is A = 1, allows 
the definition of the interaction energy density 



S^'^^x) = (-£'"*(a;)) = - 



(5i$ 



SiX{x) 



(32) 



A=l 



The single-particle potential energy density is defined as 

^P°*(a;) = lJ^d^y[^ix,y){Ti)G{y,x) + {Ti)G{x,y)i:{y,x)] (33) 

for complex fields. The local approximants to both f ™* and f p°* enter the 
energy-momentum tensor (18). 

The diagrammatic scries of $ given by Eq. (31) can be truncated at 
any level. Keeping the variational property (30), this defines a truncated 
self-consistent scheme for the KB equations (1). The so constructed self- 
energies lead to a coupling between the different species a, which obey 
detailed balance^. It has been shown [38], see also [15], that such a self- 
consistent scheme is exactly conserving at the expectation value level and 
thermodynaniic;ally c;onsistcnt at the same time. 

A prominent example is the particle-hole ring resummation in Fermi 
liquid theory in the limit of zero range four-fermion coupling with the 
following diagrams 



(34) 



E(x,j/) = 




Q 



n>2 



(35) 



(36) 



Here n counts the number of vertices in the diagram, the full dots de- 
note the external points X and {x,y). Note that compared to £™^{X) 
all coordinates are integrated in and therefore its diagrams attain an 
extra combinatorial factor 1 /n [39] . These diagrams illustrate various lev- 
els of approximation. Restricting <i> just to the one-point term provides 



'generalizing the special recipes for broad resonances given in ref. [48] 
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the Hartree approximation. From the two-point level on the collision term 
becomes finite, which also leads to finite damping widths of the particles. 
Terms, where $ has more than two internal vertices, give rise to memory 
contributions to the self-energies, cf. (10), due to the intermediate times 
of the internal vertices m'E{x,y). For this example the memory terms give 
rise to the famous InT term in the specific heat of of liquid ^He [49-51] 
at low temperatures T. 

It is possible to transcribe the variational rules to the local approxi- 
mation defining a local ^-functional replacing everywhere G by its local 
approximant, cf. (6), 



(37) 



Here X is an externally given parameter defining the reference point for 
the local approximation, and denotes the scaling factors of the vertices 
on the time or anti-time ordered branches of the contour. The variational 
rules (30) and (32) then transcribe to 



-iSjfe(X,p) 



= T 



(^i$ioc(X) 
5iG>'^{X,p) 



6iX- 



2 for real fields 

1 for complex fields 

^i$ioc(^) 



SiX+ 



(38) 
(39) 



A+=l 



In order to prove that the exact conserving properties indeed survive in 
the first-order gradient expansion (5), we shall use the conserving properties 
of the $ diagrams at each internal vertex together with the variational 
property (30). 

5.1. Properties of diagrams 

We return to the r.h.s. terms of the conservation laws (16). Due to 
the variational property (30), it is evident that they are given by closed 
diagrams of the same topology as those of <I>, however, with one point 
not contour integrated but placed on the time-ordered (— ) branch with 
coordinate X. To further reveal this relation between the "two-point" 
representation, as given by the r.h.s. of Eq. (16), and the retorded one-point 
representation of $ with respect to a chosen reference point X~ , we have 
to resolve the internal structure of the different diagrams contributing to 
Therefore, we first decompose $ into the terms of different diagrammatic 
topology 



i$loc(X) = ^i$S,(X) 
D 



(40) 



EXACT CONSERVATION LAWS 



13 



and discuss the features of any such term ^j^^.. We then enumerate the 
different vertices (r = l,...,nx) in each n\ denoting the number 
of vertices in ^(^^. For each such retarded vertex r~ wc define a one- 
point function $j^^(X;r^), which is obtained by integration and { — h} 
summation over all other vertices except r and by putting r on the time 
ordered (— ) branch at coordinate X. Subsequently we enumerate the few 
Green functions attached to r~ and label them as G~'^ or , depending 
on the line sense pointing towards or away from r~ , respectively. Thus, 
the retarded function ^^{X;r~) can be represented by a diagram, where 
the retarded point is explicitly pulled out. We draw this as a kind of 
"parachute" diagram which in local approximation becomes 



i$£,(X;r-) = 

r~ 

X (Ti)Gr'^ Pi ) • • • iCf ^ - (X, pi ) • • • (41) 

The Green functions attached to the external retarded point r~, form 
the "suspension c;ords" . The "canopy'' part C^^ spccificis the rest of the 
diagram which is not explicitly drawn^. Here all Green functions are taken 
in the local approximation, i.e. with coordinate X given by the reference 
point. From the variational principle (30) it is clear that $^(X; r~) can be 
interpreted in different equivalent ways depending on which line attached 
to r~ being opened 

i$g,(X;r-) = tJ -0^,G-\X,p)i:^r{X,p) (42) 

= ^/ (43) 

Here Gj (or Gj) is one of the suspension cords with the arrow pointing to- 
wards (or away) from r. No summation over 7 is implied by these relations! 
The self-energy terms S^*"^ and Y^^^^ are given by those subdiagrams of 

*For definiteness wc have drawn a diagram with 4 propagators linking to the retarded 
point r~ . All considerations, however, are independent of the coupling scheme which 
can even vary from vertex to vertex. For simple <E> diagrams, such as those with only two 
vertices at all, the canopy part reduces simply to a single point, cf. the second diagram 
in Eq. (34). 
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$i^^(r~) complementary to and G^, respectively. According to the 
variational rule (30) they contribute to the self-energy of a given species a 
as 

-iEL(X,p) = T-^^ = -Y.Y.Y.-^^"-^^^P)^^ (44) 

and similarly for leading to j.. The Kronecker symbol 5a-f projects 
on species a. The summation in r runs over all vertices in the ^^f. diagram. 
It is clear that each Green function in appears once in the count of 
and once in the count of This precisely matches with the two 

terms required in Eq. (15). Wc further note that the variation (32) of $ 
with respect to the coupling strength \{X) can be represented as 

^int(^) = -^^^gP = -E E *Sc(^;r-). (45) 

D reD 



5.2. Charge-current conservation 

The discussion above shows that any term of the local part of the charge 
leak Q\oc{X) in Eq. (16) is given by a closed diagram of the same topology 
as the )• Indeed, weighting $j^^(X;r^) with the charges of the in- 

and out-going Green functions reproduces piece by piece the local gain and 
loss terms in Eq. (15). The sum over all possible retarded vertices r in 
diagram D and the sum over all diagrams indeed exactly construct Qioc as 

^3loc(^) = / (04Ee«^«"Toc(^.P) 

= E E f E ^7 - E «7 J i^ioc(^; r-). (46) 

D reD \ 7 7 / 
^ V ' 

=0 

Since the charge sum vanishes at each vertex, Qioc identically vanishes. 

The fact that the local collision term part vanishes is indeed trivial. The 
point here is that the cancellation occurs diagram by diagram in terms of 
$^p(X;r~). This has the important consequence that the gradient terms 
exactly cancel out, too, i.e. 

E / (04e-^^a;7oc(^,P) = OQ^oc ^ 0, (47) 

since they arc generated by applying linear differential operations to the 
integrand of ^{^^{X;r~), while the charge factors are constants. There- 
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fore, we have verified that Q(X) = 0. This proves current conservation, 
where the current has the original Noether form. This proof applies to any 
conserved current of the underlying field theory, which relates to a global 
symmetry. 

5.3. Energy-Momentum Tensor 

In a similar way as above the four-momentum weighted terms become 

^ V ' 

=0 

xCf;..^^...(X;pi, . . . . . .){Ti)G^'HX,Pi) ■ ■ •iGf^-(X,pi) • ■ • . 

Here we have used the integrand form (41) of the parachute diagram 
$(X;r~), since now the weights are momentum dependent. This local 
collision-term part again drops out in the same way as above. However, 
the gradient correction to T'^ now involves momentum derivatives which 
can act on the p^'-factors in Eq. (14) through partial integrations. Indeed, 
only those momentum derivatives survive that act on any of the p'^-factors, 
since momentum derivatives acting on any of the Green functions leave the 
vanishing pre-factor in Eq. (48) untouched. 

The actual evaluation of the gradients is subtle and depends on the 
detailed topological structure of the diagram. For this purpose we use the 
diagrammatic rules (24) for the gradient terms with the double line as 
d^G{X,p), dashed lines as oc dS{p)/dpf'' . For the subsequent analysis we 
also introduce a new diagrammatic element, a genuine two-point function 
which gives the four-momentum ^''-factor "flowing" through the line 

' u ^ ' u ^ 
p"=J*Bi^^ with j|K5i?^= -J|l^^= -i- =5'^^ (49) 



The last relation results from partial integration and depends only on the 
sense of the dashed line. 

Note that the p'^ factors in the r.h.s. of Eq. (14) occur outside the 
gradient expression. Using the convolution theorem (28), we first pull the 
factors into the gradient terms in the following manner (using g^"^ = 

dp^/dp^) 

0{E-G}p'^-p'^0{G-S} 

= g'^'d^ (S • G + G • S) -h O {E • G • p'' - p'^ • G • E} . (50) 

^ V ' ^ V ' 
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Please, note the order of the terms under the O operator as they do not 
commute] The first term, E", arose from the Poisson bracket term in Eq. 
(28). In accord with Eq. (33), the four- momentum integration obviously 
provides the potential energy density term 

5 / (0^^i'(^'^) = g^^d.ei^^i.x) (51) 

of the energy-momentum tensor (18). 

Thus, we expect the second E2 term to generate the remaining interac- 
tion energy density part g^'^£'j"*(X) of O^"^. The four- momentum integra- 
tion of the E2 term again leads to a coalescence of the two external points. 
Thus the corresponding diagrams are of the parachute type 




where both reference points for the gradient expansion coalesce to r^, cf. 
Eqs. (42) and (43). Here the factor occurs in sequence at each of the 
suspension cords reflecting the sum over 7 and 7. In order to exploit the 
fact that thereby '}2-yP'^i ^ IZ^P'^f vanishes, all the different diagrams in the 
bracket in (52) have to be evaluated in the same way. Omitting all labels, 
one obtains for the first diagram in Eq. (52) 




(a) (b) 
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where C = ^5C/5iG and the four-scalar product between the double and 
dashed lines is implied in each diagram. The sequence (a) to (f) defines all 
diagrams resulting from the gradient expansion. The two diagrams (a) and 
(/?) shown in (54) are the only ones that finally survive the sum, i.e. the 
sum over "suspension cords" in Eq. (52). In detail: diagrams (a) and (b) 
specify the gradient terms arising from the space-time derivatives acting 
on Green functions within the canopy part. Using addition theorem (26) 
the dashed lines can be first linked from the bottom point to one definitely 
chosen upper suspension points f, as shown in diagram (a), and from there 
then further linked to the end points of the double line. However, the 
latter terms lead to p-derivatives entirely within the canopy, which finally 
drop out due to the vanishing sum. Diagrams (c) to (f) sequentially 
take the space-time gradients of the Green functions in the suspension 
cords converting them to a double line in each case. Also here only terms 
survive, where the momentum derivative acts on the p" factor, leading to 
digram (/?). 

One can now compile the terms in Eq. (52) provided point f is kept 

fixed. Then from [a) only a single term survives, namely that where the 
p^ factor is on the line linking to f, while from (/?) each term survives. The 
net result leads to a common ^(''''-factor, cf. (49), times a total derivative 
of the entire parachute diagram, i.e. 



D,reD 




= -ig^'^d, ^^foc{X,r-) (55) 



D,reD 



SiX- 

= -g^^d^S'^'iX) (56) 

from Eq. (45). Together with relation (51), this gives 

T-^{X) = -d^ [S^liX) - £C{X)) (57) 

for the r.h.s of conservation law (14). It is seen to be determined by 
the full divergence of the difference between the interaction energy and 
single-particle potential energy densities defined by Eqs. (32) and (33), 
now however evaluated in the local approximation, i.e. with no gradient 
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terms in the Wigner representation: 




(58) 





-+■ 



(59) 



Here the superscript R denotes the corresponding retarded function. For 

each diagram the contribution to £j"*(X) results from the correspond- 
ing terms of £f°^{X) just by scahng each term by the number of vertices 
n\ over the number of Green functions no- 



The quantum transport equations in the form originally proposed by 
Kadanoff and Baym, (5) or equivalently (10), have very pleasant generic 
features. As possible memory effects in the collision term are to be included 
only up to first-order space-time gradients, they are local in time to the 
extent that only the knowledge of the Green functions and their space-time 
and four-momentum derivatives at this time is required to determine the 
future evolution. They further preserve the retarded relations among the 
various real-time components of the Green function. 

In this paper we have shown that they also possess exact rather than 
approximate conservation laws, related to global symmetries of the sys- 
tem, if the scheme is ^-derivable. The same Noether currents and the 
same energy-momentum-tensor [15] as those for the original KB equations, 
however now in their local approximation forms, are exactly conserved for 
the complete gradient-expanded KB equation, i.e. the quantum transport 
equations (5). Thus, 



6. CONCLUDING REMARKS 



(X) = 0, 



with 



(60) 




(61) 




sum of single particle terms 



+g^'^{£i:',{X)-£C{X)), (62) 
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here written in terms of the product of phase-space occupation and spectral 
functions fa{X,p)Aa{X,p) = {^i)G~~^{X,p), are exact consequences of the 
equations of motion. In order to preserve this exact conserving property, 
two conditions have to be met. First, the original KB equations should 
be based on a ^-derivable approximation scheme that guarantees that the 
KB equations themselves arc conserving [38, 18]. The second condition 
is that the gradient expansion has to be done systematically, whereby it 
is important that no further approximations are applied that violate the 
balance between the different first-order gradient terms. 

Indeed, all gradient terms residing in the Poisson brackets and in the 
memory collision term cancel each other for the conserved currents such 
that the original Noether expression (61) remains conserved. This implies 
the compensation of drag- flow terms by all the other gradient terms (back- 
flow and memory flow), cf. the discussion given in ref. [18]. For the energy- 
momentum tensor the gradient terms result into the divergence of the dif- 
ference between interaction energy density and single-particle potential en- 
ergy density, cf. (62). Thereby fjj*(X) and £f°^iX) are obtained from the 
same ^-functional in the local approximation as the self-energies driving 
the equations of motion (5) and (9). The so obtained energy-momentum 
tensor is general and applies to any local coupling scheme. The energy 
component Qi^^ has a simple interpretation. The first term determines the 
single-particle energy, which consists of the kinetic and single-particle po- 
tential energy parts. Evidently this part by itself is not conserved^. Rather 
its potential energy part is compensated by the last term, i.e. £^°^{X), such 
that finally the total kinetic plus interaction energy survive. 

A typical example for the imbalance of gradient terms is the case, where 
one neglects the second Poisson bracket term in (10). This implies that 
drag-flow effects contained in the first Poisson bracket remain uncompen- 
sated. Also possible memory effects Cmem in the collision term should not 
be omitted. Otherwise the Poisson brackets |RcS^, iG'~+} + {iS~+, ReG^} 
remain uncompensated and, as a consequence, the conservation laws are 
again violated already in zero-order gradients. 

A less evident example is the modification of the gradient terms after 
the formal gradient expansion, as it has been suggested by Botermans and 
Malfliet [10], see also [18]. There one simplifies those self-energy terms 
that are involved in the Poisson brackets {ReS^, iG^'^} + {iS""*", RcGf^}, 
employing quasi-equilibrium relations. This modification implies deviations 
at second-order gradients only^, which is quite acceptable from the formal 



^Contrary to the constructions given in ref. [52]. 

^This freedom of choice is due to the fact that various redundant conilrinations of the 
KB equations lead to non-redundant equations after gradient approximation due to the 
asymmetric treatment of sums and differences of the KB equations and their adjoint 
ones in the gradient approximation. The so called mass-shell equation indeed agrees 



20 



J. KNOLL ET AL. 



point of view, when one considers a slow space-time dynamics. However, 
such kind of modifications violate the strict balance between the gradient 
terms and thus k^ad to approximate conservation laws though within first- 
order gradients. Although even in the BM case an exact conservation 
law can be formulated for an effective charge [20], which however only 
approximately coincides with the true (Nocthcr) one, exact conservation 
laws of energy and momentum could not be derived yet. 

The presence of exact conservations puts the Kadanoff-Baym formu- 
lation of quantum transport to the level of a generic phenomenological 
concept. It permits to define phenomenological models for the dynamical 
description of particles with broad damping widths, such as resonances, 
with built-in consistency and exact conservation laws, which for practical 
simulations of complex dynamical systems may even be applied in cases, 
where the smallness of the gradients can not always be guaranteed. This 
opens applications to the strong non-equilibrium dynamics of high-energy 
nuclear collisions. 

We considered here systems of relativistic bosons and/or non-relativistic 

particles with a local interaction. The generalization to relativistic fermions 
with local interactions is straightforward but involves extra complications 
resulting from the spinor structure of the kinetic equations. 

Non-relativistic systems with instantaneous interaction at finite spatial 
distance (two-body potentials) also possess exact conservation laws, if given 
by a ^-derivable approximation. Such systems can equivalently be de- 
scribed by a local field theory with mesons mediating the interactions. The 
appropriate limit towards a potential picture is obtained by treating these 
interactions instantaneously (non-relativistic limit), i.e. without retarda- 
tion. This amounts to reduce the corresponding meson Dyson equation to a 
Poisson equation over the sources of the meson fields and finally eliminating 
those meson fields. As the local field theory is conserving, the corresponding 
non-relativistic picture with potentials is conserving too. Note, however, 
that in this latter case the spatially local energy-momentum tensor does 
not exist, since energies and momenta are transfered at finite distances 
through the potentials and it is only possible to formulate the conservation 
of the total space-integrated energy and momentum. 

An important example of such a conserving approximation is given by 
the ring diagrams 



in the first-order gradient terms with the here considered transport Eq. (5), cf. [18], 
however they difi^er in higher orders of gradients. 




(63) 
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-iS(a;,y) = 




n=2 I I I I I I 



(64) 



which is just the finite range analog of the Fermi-Uquid example (34), now 
for the particle-particle rcsummation channel. Here n counts the number 
of interaction (dashed) lines representing the two-body potential V{xi — 
Xk)- The first terms in both expressions are the usual Hartree terms. The 
remaining sum leads to a conserving T-matrix type of approximation for 
the self-energies (64), which, e.g., provides thermodynamically consistent 
description of the nuclear matter [53-55]. In the dilute limit it expresses the 
self-energies through the vacuum scattering T-matrix [56,57,15]. In this 
limit it provides a collision term given by vacuum scattering cross sections 
and at the same time the gradient terms account for the appropriate virial 
corrections. The latter indirectly depend on the energy variations of the 
corresponding phase shifts, which give rise to delay time effects [58,15] and 
the corresponding changes of the underlying equation of state (energy- 
momentum tensor) [59,61-63]. Furthermore, the ,s-channel bosonization of 
the interactions in the particle- hole channel, cf. ref. [18], leads to the RPA- 
approximation. Further applications and considerations of the quantum 
transport equations will be discussed in a forthcoming paper. 

A case that still requires a separate treatment is that of derivative cou- 
pling of relativistic fields, as e.g. , in the case of the pion-nucleon interaction. 
The reason is that derivative couplings produce extra terms in the currents 
and the energy-momentum tensor, which require special treatment. 

Besides all this success at the one-particle expectation value level, one 
has to keep in mind that partial Dyson resummations still may violate the 
symmetries at the two-body correlator level and beyond. In particular, it 
means that the corresponding Ward-Takahashi identities are not necessar- 
ily fulfilled within the ^-derivable Dyson resummation scheme. On the 
other side the ^-derivable scheme provides the tools to construct the driv- 
ing terms and kernels of the corresponding higher order vertex equations 
(Bethe-Salpeter equations, etc.) which precisely recover the conservation 
laws at the correlator level [37,38,64,65]. So far such equations, however, 
could mostly be solved in drastically simplified cases (e.g., by RPA-type 
rcsummation). A particular challenge represents the inclusion of vector or 
gauge bosons into a self- consistent Dyson scheme beyond the mean-field 
level, i.e. at the propagator level, since partial Dyson resummations vi- 
olate the four-dimensional transvcrsality of the propagators. A practical 
way out of this difficulty has recently been suggested in refs. [23,64]. A 
further virtue of the ^-derivable scheme is that it apparently permits a 
renornialization of the non-perturbativc self-consistent self-energies with 
temperature- and density- independent counter terms [64]. 
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APPENDIX A 



A.l. CONTOUR MATRIX NOTATION 

In calculations that apply the Wigner transformations, it is necessary 
to decompose the full contour into its two branches — the time- ordered and 
anti-time- ordered branches. One then has to distinguish between the phy- 
sical space-time coordinates x, . . . and the corresponding contour coordi- 
nates a;'' which for a given x take two values x~ = {x~) and a;+ = (a;+) 
(/U e {0, 1, 2, 3}) on the two branches of the contour (see figure 1). Closed 
real-time contour integrations can then be decomposed as 



L 







/ dx- . 




'to 


J oo 




POO 


/ dx". 




ho 


J to 



dx'^ ... = I dx ...+ I dx+ . . . 



da;+..., (A.l) 



where only the time limits arc explicitly given. The extra minus sign of 
the anti-time-ordered branch can conveniently be formulated by a { — h} 
"metric" with the metric tensor in { — h} indices 

= i-ij) = ( ; -1 ) (A-2) 

which provides a proper matrix algebra for multi-point functions on the 
contour with "co"- and "contra" -contour values. Thus, for any two-point 
function F, the contour values are defined as 



F'^{x,y) 

F,'{x,y) 
Fij{x,y) 



F{x\y^), i,j €{-,+}, with 
a,kF>'i{x,y), F^(,x,?y) := F'''{x,y)aki 
'JikCFjiF^\x,y), a'^ = dik (A.3) 
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on the different branches of the contour. Here summation over repeated 
indices is implied. Then contour folding of contour two-point functions, 
e.g. in Dyson equations, simply becomes 

H{x\y^)=H'^{x,y) = j dz'^ F{x\ z^)G{z^ ,y^) 

= J dzF'j{x,z)G^''{z,y) (A.4) 

in the matrix notation. 

For any multi-point function the external point Xmax, which has the 
largest physical time, can be placed on either branch of the contour without 
changing the value, since the contour-time evolution from a;"^^, to x+^j, 
provides unity. Therefore, one-point functions have the same value on both 
sides of the contour. 

Due to the change of operator ordering, genuine multi-point functions 
are, in general, discontinuous, when two contour coordinates become iden- 
tical. In particular, two-point functions like iF{x,y) = (^TcA{x) B (y) 
become 



iF{x,y) = 



/ iF {x,y) iF +{x,y) \ 

\iF+-ix,y) iF++{x,y)J 
^ (rA{x)B{y)) T{B{y)A{x)) ^ 

^ {A{x)B{y) 



T-^A{x)B{y)) j 



(A.5) 



where T and T^^ are the usual time and anti-time ordering operators. 
Since there arc altogether only two possible orderings of the two operators, 
in fact given by the Wightman functions F and F^ , which are both 
continuous, not all four components of F are independent. Eq. (A.5) 
implies the following relations between non-equilibrium and usual retarded 
and advanced functions 



F«(a;,2/) := @{xo - y^) {F+- {x,y) - F-+ {x,y)) 
= F--{x,y)-F-+{x,y)=F+-{x,y) 



F++{x,y) 



F^{x,y) := -eiyo - xo) {F+-{x,y) - F-+{x,y)) 

= F--{x,y) - F+-{x,y) = F-+{x,y) - F++ {x,y),{k.Q) 

where 0(xo — j/o) is the step function of the time difference. The rules for 

the co-contour functions F etc. follow from Eq. (A. 3). 

Discontinuities of a two-point function may cause problems for differen- 
tiations, in particular, since they often occur simultaneously in products of 
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two or more two-point functions. The proper procedure is, first, with the 
help of Eq. (A. 5) to represent the discontinuous parts in F and F++ 
by the continuous F and F~^ times 0-functions, then to combine all 
discontinuities, e.g. with respect to XQ—yo, into a single term proportional 
to 6(a;o — yo), and finally to apply the differentiations. One can easily 
check that in the following particularly relevant cases 

dz {F{x\z)G{z, x^) - G{x\ z)F{z, x^)) , (A.7) 

- [ dz {F{x\ z)G{z, x^) + G{x\ z)F{z, x^)) , (A.8) 
IJ. Jc 



c 

_d_ 
dx 



- ^) ^ {F{x\z)G{z, y^) - G{x\ z)F{z, f)) 



(A.9) 



all discontinuities exactly cancel. Thereby, these values are independent of 
the placement of and x^ on the contour, i.e. the values are only functions 
of the physical coordinate x. 

For such two point functions complex conjugation implies 

{].F^+{x,y))* = iF^+{y,x) ^ iF-+(X,p) = real, 

(iF+-{x,y)y = iF+-{y,x) ^ iF+- (X, p) = real, 

(iF—{x,y)y = \F++{y,x) => {iF—{X,p)y=iF++{X,p), 

(F«(x,2/))* = F^(2/,x) ^ (F«(X,rt)* = F^(X,p), (A.IO) 

where the right parts specify the corresponding properties in the Wigner 
representation. Diagrammatically these rules imply the simultaneous swap- 
ping of all + vertices into — vertices and vice versa together with reversing 
the line arrow-sense of all propagator lines in the diagram. 
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